Historic,  archived  document 

Do  not  assume  content  reflects  current 
scientific  knowledge,  policies,  or  practices. 


ESTIMATION  OF  POND  EVAPORATION 
IN  THE  GEORGIA  COASTAL  PLAIN 

ARS-S-89 
March  1976 


AGRICULTURAL  RESEARCH  SERVICE  .  U.S.  DEPARTMENT  OF  AGRICULTURE 


Agricultural  Research  Service 
UNITED  STATES  DEPARTMENT  OF  AGRICULTURE 

in  cooperation  with 

Georgia  Institute  of  Technology 
Middle  South  Georgia  Soil  and  Water  Conservation  District 
University  of  Georgia  College  of  Agriculture  Experiment  Stations 


CONTENTS 

Page 

Abstract   

Introduction    1 

Procedures  for  estimating  evaporation    1 

Weather  Bureau  method    2 

Penman  model    - 

Modified  fractional-evaporation  equivalent  method    2 

Fractional-evaporation  equivalent  method    2 

Linear  regression  with  meteorological  factors    2 

Results  of  data  analysis   

Sensitivity  analysis    4 

Factor  analysis  of  pan  evaporation  and  meteorological  data    5 

Conclusions   " 

Literature  cited    ^ 


TABLES 

1.  Comparison  of  computed  and  observed  evaporation  rates    3 

2.  Monthly  relative  sensitivity  values  for  Tifton,  Ga   4 

3.  Correlation  matrix    ° 

4.  Eigenvalues  and  corresponding  eigenvectors    5 


i 


ESTIMATION  OF  POND  EVAPORATION 
IN  THE  GEORGIA  COASTAL  PLAIN 

By  Richard  H.  McCuen1  and  Loris  E.  Asmussen2 

ABSTRACT 

Various  techniques  to  estimate  pond  evaporation  are  compared  on  the  basis  of 
their  ability  to  reproduce  observed  evaporation  rates  from  a  floating  pan.  Methods 
that  use  only  one  meteorological  factor  provided  correlation  coefficients  that  were 
only  slightly  smaller  than  those  obtained  with  methods  requiring  values  for  two  to 
four  meteorological  factors.  Factor  analysis  indicated  that  radiation,  humidity, 
and  air  temperature  have  similar  effects  on  the  variation  in  floating-pan  evapora- 
tion; thus,  it  is  not  necessary  to  use  all  three  factors  for  predicting  evaporation 
rates.  Correlation  analysis,  sensitivity  analysis,  and  factor  analysis  indicated  that 
variation  in  air  speed  has  little  effect  on  variation  in  evaporation. 


INTRODUCTION 

The  mean  annual  rainfall  for  the  48-year  period 
1923-70  recorded  at  the  Coastal  Plain  Experiment 
Station,  Tifton,  Ga. ,  was  47.24  inches.  This  was  well 
distributed,  with  the  months  of  October  and 
November  showing  the  lowest  average  rainfall  of 
1.93  inches  and  July  the  highest  of  6.03  inches. 
However,  mean  monthly  rainfall  is  misleading  be- 
cause there  are  several  drought  periods  during  the 
crop-growing  season  each  year.  Storage  of  suffi- 
cient water  to  meet  irrigation  requirements  is  the 
obvious  solution;  hence,  estimates  of  pond  evapora- 
tion in  the  Georgia  Coastal  Plain  are  necessary  for 
estimation  of  pond  storage  requirements. 

The  net  evaporation  rate  is  controlled  by  the 
meteorological  conditions  of  the  overlying  air  and 
the  energy  state  of  the  evaporating  surface  layer. 
The  quantity  of  thermal  energy  stored  in  a  pond  or 
lake  can  also  be  a  controlling  factor  during  certain 
periods  of  the  year.  Evaporation  rates  from  open 
water  surfaces  have  been  estimated  from  mea- 
surements of  water  losses  from  pans  (16)3  and 
sunken  tanks  (1),  by  using  models  developed  from 
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energy  budgets  (13),  and  by  empirical  analysis  of 
data  (9).  Methods  that  require  measurements  of 
water  surface  properties  are  rarely  used  because 
water  surface  data  are  difficult  and  expensive  to 
collect.  Thus^  most  evaporation  estimation  proce- 
dures are  based  on  data  that  can  be  collected  above 
the  water  surface. 

Hydrologic  models  are  potentially  useful  for  pre- 
dicting the  output  from  a  hydrologic  process  and  for 
identifying  the  relative  importance  of  variables  that 
influence  the  hydrologic  process  being  modeled.  If 
the  only  objective  of  modeling  is  prediction,  it  may 
be  reasonable  to  structure  a  model  on  only  those 
variables  that  contribute  statistically  significant, 
physically  meaningful  information.  Numerous 
methods  involving  meteorological  factors  as  predic- 
tor variables  have  been  proposed  for  estimating 
evaporation  rates.  However,  past  investigations 
have  not  identified  the  marginal  value  of  different 
meteorological  factors  for  predicting  evaporation 
rates.  In  this  study  a  sensitivity  analysis,  a  factor 
analysis,  and  a  correlation  analysis  were  used  to 
examine  the  marginal  value  of  four  meteorological 
factors  (radiation,  air  speed,  relative  humidity,  and 
air  temperature)  for  predicting  evaporation  rates  in 
the  Georgia  Coastal  Plain. 

PROCEDURES  FOR 
ESTIMATING  EVAPORATION 

Numerous  methods  have  been  proposed  for  es- 
timating evaporation.  The  prediction  methods  ex- 
amined herein  are  classified  according  to  the  data 


required  to  estimate  evaporation  and  whether  or 
not  the  method  requires  empirical  calibration  of 
structural  coefficients.  The  following  three 
methods,  none  of  which  requires  the  estimation  of 
structural  coefficients,  were  examined:  the  U.S. 
Weather  Bureau4  method  (WB),  the  Penman 
method,  and  a  modified  fractional-evaporation 
equivalent  method  (MFEE).  Additionally,  the 
fractional- evaporation  equivalent  method  (FEE) 
and  four  linear  bivariate  regression  equations  were 
calibrated  from  observed  data  and  then  were  used 
to  predict  daily  evaporation. 

Weather  Bureau  Method 

The  WB  method  of  estimating  evaporation  was 
developed  by  Kohler,  Nordenson,  and  Fox  (7)  from 
an  extensive  analysis  of  data.  The  WB  method  re- 
quires observations  of  air  temperature  (T),  radia- 
tion (R),  humidity  (H),  and  air  speed  (V).  The 
method  was  originally  presented  in  nomograph 
form;  Lamoreux  (8)  provided  the  following  func- 
tional representation: 

£  =gi(T-2 12 )(0. 1024- 0.1066  In  e  R  )]  _  Q  QQQJ 

+  0.025  (e„  -  ea)°  88(0.37  +0.0041F) 

x  [0.025+(r+398.36)-2  (4.7988x  1010)  el-74a2-6"T+39S-36»]  -'(1) 

Here,  es  and  ea  are  the  vapor  pressures  at  dew  point 
temperature  and  air  temperature,  respectively. 

Penman  Model 

In  a  simplified  form,  an  energy  balance  equates 
the  net  radiant-heat  flux  density  (Rn)  to  the  sum  of 
the  latent  heat  of  evaporation  (E)  and  the  flux 
density  of  sensible  heat  (S), 

R„  =  E  +  S.  (2) 

From  aerodynamic  considerations  of  turbulent 
transfer,  the  latent  heat  flux  density  is  given  by  a 
form  of  Dalton's  law, 

E  =f0  (V)  (e0  -  ea),  (3) 

and  the  sensible  heat  transfer  is  given  by 

s  =  yf0  (V)  (T0  -  Ta),  (4) 

where  y  is  the  psychrometric  constant,  f0  (V)  is  a 
wind  function,  e  is  the  vapor  pressure,  T  is  the 
temperature,  and  the  subscripts  o  and  a  refer  to 
water  surface  and  air  quantities,  respectively. 
Equations  2,  3,  and  4  are  the  basis  for  the  Penman 
model  (13,  U),  given  by 

E0  =  (A  •  Rn  +  y  "■  Eao)  I  (A  +  y),  (5) 

where  E0  is  the  evaporation  estimate  and  A  is  the 
slope  of  the  saturation  vapor  pressure  curve  at  air 
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temperature  (Ta).  Eao  is  an  empirical  representa- 
tion of  the  latent  heat  flux  density  given  by 

Eao  =  0.35  (0.5  +  V/100)  (e0  -  ea),  (6) 

where  V  is  the  air  speed  in  miles  per  day.  Stanhill  et 
al.  (18)  showed  that  the  ratio  of  net  radiation  to 
incoming  solar  radiation  was  approximately  58  per- 
cent for  a  small  pond  in  Israel  at  approximately  the 
same  latitude  as  Tifton,  Ga.  Monteith  and  Szeicz 
(12)  calculated  a  net  radiation  ratio  of  0.53  for  an 
open  water  surface.  Da  vies  U)  indicated  that  net 
radiation  (Rn)  is,  on  the  average,  55  percent  of  the 
measured  radiation  (R).  Thus,  E0  in  equation  5  can 
be  approximated  by 

E0  =  [0.55  •  A  ■  R  +  y 

■  0.35  (0.5  +  y/100)(eo  -  ea)]/(A  +  y).  (7) 

Modified  Fractional-Evaporation 
Equivalent  Method 

The  MFEE  method  given  by  Linacre  (10)  relates 
the  fraction  of  insolation  (EIR)  used  for  evaporation 
with  the  mean  daily  ambient  temperature  Tc,  in 
degrees  centigrade,  and  the  fraction  relative  hu- 
midity H  by  the  following  empirical  nonlinear 
equation: 

EIR  =  0.02  H  Tc  +  0.9  -  H.  (8) 

Linacre  indicated  that  the  coefficient  0.02  agrees 
favorably  with  values  derived  empirically  for  evap- 
oration from  grass  surfaces.  Although  the  coeffi- 
cient for  water  surfaces  would  probably  be  smaller, 
a  value  for  water  was  not  available,  and  the  coeffi- 
cient 0.02  is  used  herein. 

Fractional-Evaporation  Equivalent 
Method 

The  FEE  method  relates  the  fraction  of  insolation 
used  for  evaporation  to  air  temperature  (Ta)  by  the 
linear  equation 

E/R=g+fTa,  (9) 

where  E  and  R  are  evaporation  and  radiation,  re- 
spectively, and/  and  g  are  regression  coefficients. 
The  FEE  method  has  been  used  by  Stephens  and 
Stewart  (19)  in  southern  Florida  and  Lane  (9)  in  the 
Western  United  States.  Values  of  the  regression 
coefficients  were  determined  empirically  for  this 
study  from  the  analysis  of  radiation,  temperature, 
and  pan  evaporation  data. 

Linear  Regression  With  Meteorological 
Factors 

The  following  linear  regression  equation  was 
used  to  relate  individually  each  of  the  four 
meteorological  factors  (radiation,  air  speed,  air 
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temperature,  and  relative  humidity)  to  daily  pan 
evaporation  rate  E: 

E=a+bF,  (10) 

where  F  is  the  meteorological  factor  and  a  and  b  are 
the  regression  coefficients.  The  regression  coeffi- 
cients were  derived  from  the  observed  evaporation 
from  a  floating  pan  as  the  criterion  variable. 

RESULTS  OF  DATA  ANALYSIS 

The  Southeast  Watershed  Research  Unit  (Ag- 
ricultural Research  Service)  is  conducting  hydro- 
logic  studies  at  a  pond  on  the  J.  B.  Walker  farm,  5 
miles  southwest  of  Tifton,  Ga.  Walker  Pond  has  a 
maximum  storage  of  9.81  acre-feet  and  a  maximum 
surface  area  of  2.56  acres.  Daily  evaporation  is 
measured  from  a  pan  floating  on  the  pond.  The 
floating  pan  is  partially  shielded  by  a  group  of  trees. 
The  average  tree  height  in  the  woodland  (pine)  east 
|  of  the  evaporation  pan  is  25  feet.  There  are  about 
650  trees  per  acre  in  this  plot. 

Meteorological  data,  including  relative  humidity, 
;  air  temperature,  air  speed,  and  incoming  radiation, 
are  measured  at  the  U.S.  Weather  Service  office  in 
Tifton,  Ga.  During  fall  and  winter,  the  prevailing 
i  winds  are  from  the  northwest  and  northeast.  Dur- 
ing the  remainder  of  the  year,  the  prevailing  winds 
are  from  the  west  and  southwest.  During  a 
33-month  period,  392  daily  measurements  of  evap- 
oration from  the  floating  pan  and  the  four 
meteorological  factors  from  the  Weather  Service 
office  were  made. 

Estimates  of  daily  evaporation  were  derived  by 
each  of  the  previously  described  methods  from  the 
392  daily  measurements  of  the  4  meteorological  fac- 
tors and  floating-pan  evaporation  data.  The  com- 
puted evaporation  rates  were  compared  with  the 
|  corresponding  measured  values.  The  Pearson 
product-moment  correlation  coefficient  R  was  used 
as  a  measure  of  the  strength  of  the  linear  association 
between  the  observed  and  computed  evaporation 
'  rates.  The  coefficient  of  determination  R2  is  the 
proportion  of  the  variance  in  the  observed  evapora- 
I  tion  rates  that  can  be  attributed  to  its  association 
with  the  computed  rates.  The  resulting  values  of  R 
and  R2  are  given  in  table  1. 

The  bivariate  linear  regression  equations  involv- 
ing humidity,  radiation,  air  temperature,  and  air 
speed  explained  38.5,  32.2,  31.2  and  1.5  percent  of 
the  total  variation  in  the  observed  daily  evapora- 
tion, respectively.  The  inadequacy  of  the  air  speed 
measurements  for  estimating  evaporation  is  dem- 
onstrated by  the  negative  value  of  the  correlation 


TABLE  1. — Comparison  of  computed  and  observed 
evaporation  rates 


Procedure  «  " 

Weather  Bureau   0.650  0.423 

Penman   626  393 

Modified  fractional-evaporation  equivalent  623  .388 

Fractional-evaporation  equivalent   579  .334 

Linear  regression:  Predictor  variable: 

Humidity   --620  .385 

Radiation   567  322 

Temperature   559  .312 

Air  speed  —  .122  .015 


coefficient.  A  negative  value  of  R  for  such  a  rela- 
tionship indicates  that  increases  in  air  speed  are 
accompanied  by  decreases  in  evaporation.  Such  a 
result  would  not  be  expected. 

The  FEE  method,  which  uses  two  meteorological 
factors  (radiation  and  air  temperature)  to  estimate 
evaporation,  explained  33.4  percent  of  the  total 
variation,  only  1.2  and  2.2  percent  greater  than 
the  percentage  of  variation  explained  by  the  bi- 
variate regression  equations  for  radiation  and  air 
temperature.  This  suggests  that  radiation  and  air 
temperature  are  intercorrelated  and  tends  to  ex- 
plain similar  variation  in  evaporation  rates. 

The  MFEE  method,  equation  8,  explained  38.8 
percent  of  the  total  variation  in  the  observed  evap- 
oration rates,  which  is  only  0.3  percent  greater  than 
the  value  resulting  from  the  bivariate  regression 
equation  for  humidity.  However,  the  coefficients 
for  the  bivariate  regression  equation  involving 
humidity  were  derived  from  data  collected  at  Tif- 
ton, Ga.,  while  the  coefficients  for  the  MFEE 
method  were  from  the  developments  of  a  previous 
investigator  (10).  Thus,  the  difference  of  0.3  percent 
in  explained  variation  is  not  necessarily  an  indica- 
tion that  humidity,  air  temperature,  and  radiation 
explain  similar  variation  in  evaporation.  Also,  the 
coefficients  of  the  MFEE  method  are  possibly  in- 
adequate for  estimating  evaporation  from  ponds  in 
the  southern  Coastal  Plain. 

The  WB  and  Penman  methods,  which  require 
values  of  four  meteorological  factors,  explained  42.3 
and  39.3  percent  of  the  total  variation,  respectively. 
The  MFEE  method,  which  involves  only  three 
meteorological  variables  and  does  not  require  cali- 
bration of  structural  coefficients,  explained  38.8 
percent  of  the  variation.  Thus,  air  speed,  which  is 
used  with  the  WB  and  Penman  methods  but  not 
the  MFEE  method,  does  not  significantly  increase 
the  percentage  of  explained  variation.  Such  a  re- 
sult was  expected,  since  the  correlation  coefficient 
for  the  bivariate  regression  equation  involving 
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air  speed  was  very  small  and  not  statistically 
significant. 

Even  though  the  structural  coefficients  for  the 
bivariate  regression  equation  involving  humidity 
were  derived  from  observed  data,  this  equation  ex- 
plained only  3.8  and  0.8  percent  less  variation  than 
the  WB  and  Penman  methods,  respectively.  Thus, 
it  is  reasonable  to  examine  the  possibility  that 
measurements  of  radiation,  air  temperature,  and 
humidity  explain  similar  variation  in  evaporation 
measurements  and  that  the  collection  of  data  for 
all  three  meteorological  factors  represents  a  dupli- 
cation of  information  and  effort  for  prediction  of 
evaporation  rates;  a  factor  analysis  will  be  used 
to  examine  the  intercorrelations  of  the  meteoro- 
logical variables  in  a  later  section. 

SENSITIVITY  ANALYSIS 

Correlation  coefficients  are  used  to  measure  the 
degree  of  association  between  variables.  The  Pear- 
son product-moment  correlation  coefficient  is  de- 
fined by 

1=1       i=i  1=1 

where  x=X—X,  y  =  Y—Y,  andn  is  the  number  ofx,y 
pairs.  The  square  of  the  correlation  coefficient  rep- 
resents the  proportion  of  the  variance  in  a  criterion 
Y  that  can  be  attributed  to  its  linear  regression  on  a 
predictor  X. 

Regression  analysis  is  used  to  define  the  relation- 
ship between  variables.  There  is  a  strong  structural 
similarity  between  correlation  analysis  and  regres- 
sion analysis  (1 7).  The  linear  regression  coefficient  b 
can  be  computed  from  the  correlation  coefficient 
and  the  standard  deviation  of  X  and  Y: 

b=  £  xiVi  I  %x*=RSJSx.  (12) 

i=l  i=i 

Equation  12  suggests  that  for  values  of  R  near  1, 
changes  in  X  will  produce  comparatively  large 
changes  in  Y. 

The  sensitivity  of  evaporation  to  a  meteorological 
factor  is  defined  as  the  rate  of  change  in  evaporation 
with  respect  to  change  in  the  meteorological  factor. 
For  a  linear  regression  equation  relating  evapora- 
tion E  and  a  meteorological  factor  F,  the  sensitivity 
of  E  to  variation  in  F  can  be  determined  by  taking 
the  derivative  of  E  with  respect  to  F: 

dEldF=b.  (13) 

Equations  12  and  13  show  that  there  is  a  direct 


relationship  between  linear  regression  and  sensitiv- 
ity. Therefore,  estimates  of  sensitivity  can  be  used 
to  determine  a  priori  the  factors  that  might  account 
for  the  most  variation  in  the  dependent  variables. 

The  value  of  sensitivity  computed  by  the  deriva- 
tive approach  of  equation  13  is  not  invariant  to  the 
magnitude  of  either  the  dependent  or  independent 
variable.  Relative  sensitivity  Rs  is  defined  as  the 
ratio  of  the  percentage  change  in  evaporation 
(AE/E)  to  the  percentage  change  in  the  magnitude 
of  the  meteorological  factor  (AF/F): 

R,=  *E/E  =  (&E\(F\  (14) 
AF/F  \AF/\E/ 

Since  the  relative  sensitivity  is  a  function  of  F 
and  E ,  it  can  be  used  for  comparing  the  sensitivity 
of  the  different  meteorological  factors  that  in- 
fluence evaporation.  A  negative  value  of  Rs 
indicates  that  as  the  value  of  the  meteorological 
factor  increases  evaporation  will  decrease. 

In  order  to  relate,  relative  sensitivity  values  to  the 
influence  of  the  meteorological  factors  in  estimating 
evaporation  rates,  it  is  necessary  to  specify  a  model 
of  the  evaporation  process  from  which  Rs  values  can 
be  derived.  Use  of  a  model  based  on  the  underlying 
physical  principles  will  provide  Rx  values  that  are 
indicative  of  the  relative  influence  of  variables  in  the 
evaporation  process.  Thus,  the  Penman  model  was 
selected  because  it  is  based  on  a  simplified  energy 
balance  equation.  With  average  monthly  values  (20) 
of  meteorological  factors  and  the  Penman  model, 
monthly  values  of  relative  sensitivity  were  deter- 
mined (table  2). 

The  mean  annual  relative  sensitivity  values  for 
the  four  meteorological  factors  indicate  that  humid- 

TABLE  2. — Monthly  relative  sensitivity  values  for 
Tifton,  Ga. 


Relative  sensitivity 
Mean  daily  of  the  meteorological  factor 


Month 

evaporation 
(inch/day) 

Humidity 

Aii- 
speed 

Radiation 

Temperature 

Jan.    .  . 

0  082 

-1.034 

0.256 

0.655 

0.661 

Feb. 

101 

-.811 

.258 

.652 

.688 

Mar 

132 

-.672 

.224 

.698 

.637 

Apr. 

172 

-.569 

.199 

.732 

.621 

May 

196 

-.542 

.189 

.745 

.627 

June 

199 

-.521 

.158 

.777 

.589 

July 

184 

-.641 

.128 

.819 

.484 

Aug.  .. 

180 

-.656 

.131 

.815 

.495 

Sept. 

148 

-.763 

.152 

.785 

.554 

Oct. 

119 

-.805 

.189 

.732 

.629 

Nov.  ,  . 

092 

-.846 

.232 

.671 

.677 

Dec. 

...  .073 

-1.184 

.277 

.626 

.717 

Mean 

-.753 

.199 

.726 

.615 
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jity  is  the  most  sensitive  factor  and  that  radiation 
and  temperature  are  slightly  less  sensitive.  The 
!  mean  Rs  value  for  air  speed  indicates  that  a 
.  10-percent  change  in  air  speed  produces  a  2-percent 
i  change  in  evaporation.  Linsley,  Kohler,  and 
,Paulhus  (11)  indicated  that  a  10-percent  change  in 
air  speed  causes  a  change  in  evaporation  of  1  to  3 
percent.  Thus,  the  relative  sensitivity  values  de- 
i  rived  herein  appear  reasonable.  A  ranking  of  the 
jmean  sensitivity  of  the  four  factors  corresponds 
f  exactly  to  the  ranking  of  correlation  coefficients  for 
the  bivariate  regression  equations  of  table  1.  In 
j  both  rankings,  there  is  a  very  noticeable  difference 
|  between  air  speed  and  the  three  remaining  me- 
teorological factors. 

The  monthly  Rs  values  of  table  2  indicate  that 
evaporation  rates  are  most  sensitive  to  variation  in 
humidity  during  winter.  During  summer,  variation 
in  radiation  produces  a  greater  change  in  evapora- 
tion than  does  variation  in  humidity.  From  April  to 
June,  both  temperature  and  radiation  are  more  sen- 
sitive than  humidity.  Thus,  the  values  of  relative 
sensitivity  are  potentially  useful  for  indicating  the 
relative  importance  of  the  factors  that  control  the 
evaporation  process  and  the  contribution  of  each 
i  factor  in  the  reproduction  of  measured  evaporation 
rates. 

FACTOR  ANALYSIS  OF  PAN 

EVAPORATION  AND 
METEOROLOGICAL  DATA 

Factor  analysis,  a  multivariate  statistical  tech- 
nique, can  be  used  to  examine  the  intercorrelations 
between  the  observed  meteorological  and  pan 
:  evaporation  data.  Such  an  analysis  is  desirable  be- 
i  cause  it  can  lead  to  a  reduction  in  the  number  of 
"independent"  variables  required  to  predict  a  "de- 
pendent" variable.  Elimination  of  predictor  vari- 
ables that  explain  similar  variation  in  the  criterion 
can  reduce  the  cost  of  future  data  collection.  A  re- 
duction in  the  number  of  predictor  variables  may 
reduce  prediction  ability  somewhat,  but  for  vari- 
ables that  are  highly  intercorrelated,  the  reduction 
is  often  insignificant. 

Factor  analysis  transforms  the  correlation  matrix 
of  the  criterion  and  predictor  variables  into  a  set  of 
variates  called  factors.  These  factors  are  linear 
combinations  of  the  original  variables  and  have  the 
advantage  of  being  uncorrelated.  The  coefficients  of 
the  factors  are  called  factor  loadings. 

The  correlation  matrix  of  table  3  was  derived 
from  392  daily  measurements  of  four  meteorological 
factors  and  evaporation  from  a  floating  pan.  With 


the  exception  of  associations  involving  air  speed, 
the  correlation  coefficients  between  each  pair  of 
variables  exceeded  0.5.  The  correlation  coefficients 
between  air  speed  and  each  of  the  other  four  vari- 
ables are  between  zero  and  -0.3.  But  the  correla- 
tion matrix  by  itself  does  not  provide  the  means  for 
determining  whether  or  not  the  four  variates  other 
than  air  speed  explain  similar  information. 

A  factor  analysis  of  the  correlation  matrix  of  table 
3  produced  the  eigenvalues  and  factor  pattern  of 
table  4.  In  the  first  factor  all  of  the  factor  loadings 
are  large  except  the  factor  loading  for  air  speed.  The 
factor  loading  for  air  speed  is  large  in  the  second 
factor  while  the  other  four  factor  loadings  are  small. 
In  the  remaining  three  factors,  all  the  loadings  are 
comparatively  small. 

Table  3. — Correlation  matrix 

Pan  Air  Radia- 
Variable          evap.    Temp.      speed       tion  VPD 
 (1)        (2)  (3)  (4)  (5) 

(1)  Pan 

evaporation  ..  1.0000    0.5587    -  0.1224    0.5667  0.6203 


(2)  Temperature .. .  .5587  1.0000  -.2914  .5282  .8240 

(3)  Air  speed  -.1224  -.2194  1.0000  -.1981  -.2939 

(4)  Radiation  5667  .5282  -.1981  1.0000  .7237 

(5)  Vapor  pressure 

deficit 

(VPD) .  6203  .8240  -.2939  .7237  1.0000 


Table  4. — Eigenvalues  and  corresponding 
eigenvectors 

Eigenvalues  for  eigenvector  of — 

Variable  3.01  093  ^049  !044  M3 
 (1)       (2)        (3)         (4)  (5) 

(1)  Pan  evaporation  .. .  0.779  0.224  -0.305  -0.491  0.002 

(2)  Temperature  854  .057  .478  -.051  .190 

(3)  Air  speed  -.358  .925  .089  .085  -.021 

(4)  Radiation  817  .101  -.366  .420  .110 

(5)  Vapor  pressure 

deficit  937     .011       .171       .121  -.279 


Although  methods  exist  for  determining  the 
number  of  statistically  significant  factors  (2,  3,  15), 
sl  universally  accepted  standard  for  when  to  stop 
factoring  has  not  been  developed  (5).  Furthermore, 
information,  even  though  statistically  significant, 
may  have  no  practical  value  (5).  Therefore,  a 
method  proposed  by  Kaiser  (6)  through  extensive 
data  analysis  was  used  to  determine  the  number  of 
significant  factors.  Kaiser  recommended  that  the 
factors  corresponding  to  eigenvalues  of  the  correla- 
tion matrix  that  are  greater  than  1  could  be  consid- 
ered significant. 

The  eigenvalues  of  table  4,  determined  by 
Kaiser's  method,  indicate  that  only  the  first  factor  is 
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significant.  The  large  factor  loadings  in  the  first 
factor,  which  correspond  to  radiation,  vapor  pres- 
sure deficit,  air  temperature,  and  pan  evaporation, 
suggest  that  observations  of  the  three  meteorologi- 
cal and  pan  evaporation  variables  contain  similar 
information  even  though  the  four  variables  are  dis- 
similar. If  the  three  meteorological  factors  explain 
similar  variation  in  pan  evaporation,  then  a  predic- 
tion equation  including  more  than  one  of  the  three 
meteorological  factors  should  not  be  expected  to 
provide  significantly  greater  correlation  than  a  pre- 
diction equation  involving  only  one  of  the  variables. 
Comparison  of  the  correlation  coefficients  of  table  1 
confirms  the  findings  of  the  factor  analysis  and  indi- 
cates that  the  methods  involving  more  than  one 
meteorological  factor  provide  correlation  coeffi- 
cients only  slightly  higher  than  those  obtained  with 
the  bivariate  regression  equations.  Furthermore, 
the  very  low  correlation  between  air  speed  and 
evaporation  supports  Kaiser's  empirical  rule  for  de- 
termining the  number  of  significant  factors.  The 
sensitivity  analysis  also  indicated  that  variation  in 
air  speed  had  little  effect  on  variation  in  evaporation 
rates. 


CONCLUSIONS 


Since  true  rates  of  evaporation  from  the  pond 
surface  were  not  available,  the  observed  floating- 
pan  evaporation  rates  are  used  to  represent  them. 
Such  an  assumption  seemed  reasonable  since  the 
mean  annual  evaporation  from  the  floating  pan  was 
only  about  9  percent  less  than  the  annual  rate  re- 
ported for  the  Tifton  area  by  the  U.S.  Department 
of  Commerce  (20).  During  fall  and  winter,  air  moves 
primarily  out  of  the  northwest  and  northeast.  The 
wooded  area  north  of  Walker  Pond  could  signifi- 
cantly influence  evaporation  rates  and  thus  account 
for,  in  part,  the  difference  between  the  observed 
and  computed  mean  annual  evaporation.  The  mean 
observed  May-to-October  evaporation  expressed  as 
a  percentage  of  the  total  annual  evaporation  dif- 
fered from  the  reported  evaporation  by  only  2 
percent.  From  May  to  October  the  wind  moves 
primarily  from  the  west  and  southwest.  During 
this  period  the  wooded  area  on  the  northern  bank 
of  the  pond  would  not  significantly  influence  evap- 
oration rates.  Furthermore,  a  recent  study  (1)  in- 
dicated that  the  average  ratio  of  evaporation  from 
Lake  Hefner  to  evaporation  from  tanks  with  diam- 
eters of  9  and  15  feet  and  located  approximately 
100  feet  from  the  lake  is  1.01. 

As  indicated  previously,  an  understanding  of  the 
evaporation  process  can  provide  the  means  for  ob- 
taining better  estimates  of  evaporation  rates.  Such 
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estimates  of  evaporation  can  then  be  used  to  pro- 
vide more  accurate  estimates  of  the  storage  re- 
quired for  irrigation  supplies  during  drought.  A 
sensitivity  analysis  was  used  to  examine  the  rela- 
tive importance  in  the  Georgia  Coastal  Plain  of 
four  meteorological  variables  thought  to  influence 
evaporation.  Mean  annual  relative  sensitivity 
values  indicated  that  humidity,  radiation,  and  air 
temperature  are  of  similar  importance,  while  vari- 
ation in  air  speed  has  a  much  smaller  influence  on 
evaporation  rates.  However,  mean  monthly  rela- 
tive sensitivity  values  indicated  that  humidity  is 
the  controlling  factor  in  winter,  while  radiation 
exerts  a  greater  influence  on  evaporation  rates  in 
summer. 

The  correlation  analysis  of  table  1  indicates  that 
the  WB  method  is  potentially  capable  of  providing 
the  best  estimates  of  evaporation.  If  a  method  for 
estimating  net  radiation  is  available,  the  Penman 
method  also  provides  reasonable  estimates  of  daily 
evaporation.  However,  the  correlation  coefficients 
obtained  with  three  of  the  bivariate  regression 
equations  were  only  slightly  less  than  the  correla- 
tion coefficients  obtained  with  the  WB  and  Penman 
equations;  this  suggests  that  the  use  of  more  than 
one  meteorological  factor  is  redundant  when  the 
objective  of  modeling  is  prediction. 

Whereas  the  coefficients  of  the  WB  and  Penman 
methods  were  derived  at  other  locations,  the  re- 
gression coefficients  for  the  bivariate  regression  , 
equations  were  derived  from  data  collected  at  Tif- 
ton, Ga.  Therefore,  a  factor  analysis  was  used  to 
examine  the  intercorrelations  between  the  me- 
teorological factors  and  floating-pan  evaporation 
rates.  The  factor  analysis  supported  the  conclusions 
of  the  correlation  analysis.  Specifically,  the  factor 
analysis  indicated  that  radiation,  air  temperature, 
vapor  pressure  deficit,  and  pan  evaporation  explain 
similar  variation,  and  thus,  for  the  Georgia  Coastal 
Plain,  the  use  of  more  than  one  of  the  three  me- 
teorological factors  (vapor  pressure  deficit,  radi- 
ation, and  air  temperature)  may  be  redundant  for 
predicting  evaporation. 

The  correlation  analysis,  sensitivity  analysis,  and 
factor  analysis  indicated  that  evaporation  rates  are 
only  slightly  influenced  by  variation  in  air  speed. 
The  sensitivity  analysis  showed  that  a  variation  in 
air  speed  of  10  percent  results  in  a  2-percent  change 
in  evaporation.  Thus,  the  very  low  correlation  coef- 
ficient, -0. 122,  between  airspeed  and  pan  evapora- 
tion should  have  been  expected.  The  factor  analysis 
also  indicated  that  air  speed  was  not  a  significant 
factor  in  estimating  floating-pan  evaporation. 


LITERATURE  CITED 

(1)  Allen,  J.  B.,  and  Crow,  F.  R.  1971.  Predicting  lake  evapora- 

tion by  performance  of  evaporation  ponds,  pans  and 
tanks.  Trans.  ASAE  14(3):  458-463. 

(2)  Bartlett,  M.  S.  1950.  Tests  of  significance  in  factor  analysis. 

Br.  J.  Psychol. ,  Stat.  Sect.  3:  77-85. 

(3)   1951.  A  further  note  of  tests  of  significance  in  factor 

analysis.  Br.  J.  Psychol.,  Stat.  Sect.  4: 1-2. 

(4)  Davies,  J.  A.  1967.  A  note  on  the  relationship  between  net 

radiation  and  solar  radiation.  Q.  J.  R.  Meteorol.  Soc.  93: 
109-116. 

(5)  Harman,  H.  H.  1960.  Modern  factor  analysis.  474  pp.  Uni- 

versity of  Chicago  Press,  Chicago. 

(6)  Kaiser,  H.  F.  1960.  The  application  of  electronic  computers 

to  factor  analysis.  Educ.  Psychol.  Meas.  20(1):  141-151. 

(7)  Kohler,  M.  H.,  Nordenson,  T.  J.,  and  Fox,  W.  E.  1955. 

Evaporation  from  pans  and  lakes.  U.S.  Weather  Bur. 
Res.  Pap.  38,  53  pp. 

(8)  Lamoreux,  W.  M.  1962.  Modern  evaporation  formulas 

adapted  to  computer  use.  Mon.  Weather  Rev.,  vol.  90, 
January. 

(9)  Lane,  R.  K.  1964.  Estimating  evaporation  from  insolation. 

J.  Hydraul.  Div.,  Proc.  Am.  Soc.  Civ.  Eng.  90(HY5): 
33-11. 

(10)  Linacre,  E.  T.  1967.  Climate  and  the  evaporation  from 
crops.  J.  Irrig.  Drain.  Div.,  Proc.  Am.  Soc.  Civ.  Eng.  93 
(IR4):  61-79. 


(11)  Linsley,  R.  K.,  Kohler,  M.  A.,  and  Paulhus,  J.  L.  H.  1958. 

Hydrology  for  engineers.  340  pp.  McGraw-Hill,  New 
York. 

(12)  Monteith,  J.  L. ,  and  Szeicz,  G.  1962.  Radiative  temperature 

in  the  heat  balance  of  natural  surfaces.  Q.  J.  R.  Meteorol* 
Soc.  88:  496-507. 

(13)  Penman,  H.  L.  1948.  Natural  evaporation  for  open  water, 

bare  soil,  and  grass.  Proc.  R.  Soc.  London,  Ser.  A,  193: 
120-145. 

(14)   1956.  Estimating  evaporation.  Trans.  Am.  Geophys. 

Union  37:  43-46. 

(15)  Rao,  C.  R.  1955.  Estimation  and  tests  of  significance  in 

factor  analysis.  Psychology  20:  93-111. 

(16)  Shulyakovskiy,  L.  G.  1969.  Formula  for  computing  evap- 

oration with  allowance  for  the  temperature  of  the  free 
water  surface.  Sov.  Hydrol.  Sel.  Pap.  6:  566-573. 

(17)  Snyder,  W.  M.  1971.  Summary  and  evaluation  of  methods 

for  detecting  hydrologic  changes.  J.  Hydrol.  12: 311-338. 

(18)  Stanhill,  G.,  Hofstede,  G.  W.,  and  Kalma,  J.  D.  1966.  Radi- 

ation balance  of  natural  and  agricultural  vegetation.  Q.  J. 
R.  Meteor.  Soc.  92:  128-140. 

(19)  Stephens,  J.  C,  and  Stewart,  E.  H.  1963.  A  comparison  of 

procedures  for  computing  evaporation  and  evapotranspi- 
ration.  Int.  Assoc.  Sci.  Hydrol.  Pub.  62,  pp.  123-133. 

(20)  U.S.  Department  of  Commerce.  1968.  Climatic  atlas  of  the 

United  States.  104  pp.  Washington,  D.C. 


*1796-G. P. 0.-1750-S/671 -583/46 


7 


.  S.  DEPARTMENT  OF  AGRICULTURE 

AGRICULTURAL  RESEARCH  SERVICE 
SOUTHERN  REGION 
P.  O.  BOX  5  3326 
NEW  ORLEANS,  LOUISIANA  70153 


OFFICIAL  BUSINESS 

PENALTY  FOR  PRIVATE  USE,  S300 


POSTAGE  AND  FEES  PAID 
U.  S.  DEPARTMENT  OF 
AGRICULTURE 
AGR  101 


